Statistical Properties of The First Excited State of an Interacting Many Particle 

Disordered System 



Richard Berkovits, 1 Yuval Gefen, 2 ' 1 Igor V. Lerner, 3 and Boris L. Altshuler 4,5 

1 The Minerva Center, Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel 
Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 
'School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, United Kingdom 
4 NEC Research Institute, 4 Independence Way, Princeton, New Jersey 08540 
5 Department of Physics, Princeton University, Princeton, New Jersey 08545 
(Dated: Mar. 3, 2003, version 3.1) 

The distribution Pi of the first many-body excitation energy of a weakly and moderately inter- 
acting electron gas in a finite conductor (in the diffusive regime) is calculated. As the interaction is 
increased, Pi crosses over from Wigner-Dyson to Poisson. We characterize this transition through 
the inverse participation ratio in Hilbert space, and examine its manifestation in a projected 2- 
dimensional space. 

PACS numbers: PACS numbers: 71.30.+h, 64.60.Ak, 73.20.Fz 



One of the most important physical quantities used 
to characterize finite conductors is energy level statistics. 
Not only does it provide an efficient mathematical frame- 
work to study such systems, but it is also correlated with 
other physical observables (conductance, spectral corre- 
lations, absorption spectrum, orbital magnetization, heat 
capacity when the number of electrons is fixed, etc.). 
This, in principle, allows for the "measurement of level 
statistics". In weakly disordered conductors ("diffusive 
disorder", as well as in chaotic cavities), it is known that 
the single electron energy levels obey the Wigner-Dyson 
(WD) statistics. 

One can further investigate the many-body spectra of 
such systems. Let us first take the interaction strength 
to vanish. In this case the distribution of the normal- 
ized energy gap between the ground-state (|0)) and the 
first-excited (|1)), s = (Ex - E )/(E 1 - E ), follows the 
WD distribution (a single particle-hole pair at the Fermi 
energy is involved). Hereafter (...) denotes an average 
taken over disorder realizations. Once we consider the 
gap between the m th and the to + 1 st many-body states 
with to » 1, we expect the statistics to follow the Pois- 
son distribution. This is due to the fact that, typically, 
the m th and the m+ 1 st many-body states correspond to 
mutually very different particle and hole configurations. 
Since they are not connected with each other through 
a single-particle operator, level repulsion is greatly re- 
duced. Once we turn on electron-electron interaction, we 
expect the gap between the m th and the m + 1 st many- 
body states (to >> 1) of such a system to satisfy the 
WD distribution. Indeed, the motivation behind the in- 
troduction of random matrix theories was to account for 
high-lying spectra of interacting systems. Table 1 sum- 
marizes the existing knowledge regarding the statistics of 
the excitation spectrum of a non-integrable electron gas. 
For the sake of definiteness we shall consider hereafter 
spinless electrons moving under the influence of weak 
disorder which renders the motion of a single electron 
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TABLE I: The many-body spacing statistics of diffusively sys- 
tems: the item marked by ? is found to cross-over ( with 
increasing interaction strength) from WD to Poisson 



diffusive. 

The item in Table 1 which is marked by a question 
mark represents the distribution of the first excitation 
energy of a (weakly or moderately weakly) interacting 
electron gas, -Pi(s). This quantity is indeed the main 
object of our investigation. Earlier studies |]j of P% have 
addressed different parameter regimes, or were at times 
inconclusive. 

The main results of our analysis are: 

1. We find numerically that as the electron-electron 
interaction strength is increased, the statistics of 
the first excitation energy crosses over from WD to 
Poisson (cf. Fig. QJ. We characterize this transi- 
tion quantitatively (cf. Fig[21and Eq^J. 

2. We relate the WD-to-Poisson transition to the 
statistics of the diagonal and off-diagonal elements 
of an effective 2x2 matrix, depicting the mixing 
of the many body-sates |0) and |1) (Fig|S]and Eqs. 

EEEJ. 

3. We show how the WD-to-Poisson crossover is man- 
ifest in the distribution of the participation ratio, 
i.e., the structure of the many-body wave function 
in Hilbert space. 

4. Finally, and very importantly, we show that as the 
interaction is increased the first excited state, which 
originally (at zero interaction) corresponds to a sin- 
gle particle-hole excitation, involves an increasing 
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number of particle-holes. But this is up to a cer- 
tain value of the interaction (r s ss 1). Upon fur- 
ther increase of the interaction strength this trend 
reverses. This picture is obtained by computing 
(numerically) the inverse participation ratio of the 
ground and first-excited states (in Hilbert space) as 
function of the interaction (cf. Fig. |3J) 

The model. We consider the following tight-binding 
Hamiltonian, describing spinless interacting electrons 
moving in the presence of a random potential: 

# = E 6k ^ nk 'i - v Yl\ a h+i a k,3 + a l+i,j a k,j + h.c] 
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The exact many-body states may be expressed in terms 
of Fock determinants : 
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Here {k,j} denotes a lattice site, the number operator 
is rifej and e k j is the site energy, selected randomly and 
uniformly over the interval [— W/2, W/2]. We study two- 
dimensional arrays with periodic boundary conditions. 
Only nearest-neighbor interactions U are considered 2, 
3]. " 

We present numerical results pertaining to N = 4 — 8 
electrons residing on 4 x 4 lattices with M = 16 sites 
(u = N/M). Data for larger systems appeared to 
be compatible with our results. The disorder strength 
W = 8V was chosen so that the single-particle localiza- 
tion length is larger than the system size, with spectral 
correlations following the Gaussian orthogonal ensemble 
(GOE) on small energy scales (i.e., for energy differences 
Ae < h/t flight ; t flight is the time-of- flight of the electron 
throughout the system). The motion of the electron is 
diffusive or marginally diffusive, l e i < L, with L being the 
quantum dot's linear size and l e i - the elastic mean- free- 
path.) Varying the electron density in the experiment 
amounts to varying the dimensionless parameter r s , 
We have carried out exact diagonalization of the many- 
particle Hamiltonian using the Lanczos method, and ob- 
tained the eigenvalues E a and eigenvectors \a), a = 0, 1. 

To establish a useful frame of reference for discussing 
the results of our analysis, it is convenient to study the 
self-consistent Hartree-Fock (SCHF) Hamiltonian corre- 
sponding to this system. This Hamiltonian reads 

Hhf = E e k,jn ktJ - V^2[al d+1 a k ,j + at+i,j a k,j + h.c] 

k,j k,j 

~J2nk,jU{n k ±i t j±i)o -^2a f k:j a k±1 j ±1 U(al ±1:j±1 a kj )o,(2) 

k,j k,j 

Here (. . .)o denotes a quantum average taken over the 
ground state, the latter being calculated self-consistently. 
The SCHF single-electron eigenvectors {IV 1 ™)} and eigen- 
values {ei] are obtained through a self-consistent diago- 
nalization of the Hartree-Fock Hamiltonian. 



where the sum is over all possible permutations of N 
states out of M with M > In > in-i > ■ ■ ■ > h > 
h > lj cjj is the creation operator of the n th single- 
particle HF state and |vac) is the vacuum state. For the 
non-interacting case, as well as for the case for which 
the SCHF provides an exact solution, the ground state 
(the first excited state) corresponds to C\^,...,N — 1 
(Ci,2,...,at-i,at+i = 1) while all other coefficients are 
equal to zero. This is compatible with a Koopmans-like 
picture. In this case the participation ratio (see below) 
is equal to 1. 

An important quantity which is particularly suitable to 
describe scenarios intermediate between the WD (77 = 0) 
and the Poisson (77 = 1) limits is p 
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WD to Poisson crossover. The distribution of s as 
function of the interaction strength for 1000 different re- 
alizations of disorder is plotted in Fig. ^ It is seen that 
Pi(s) crosses over from a Poisson-like to a WD-like dis- 
tribution as the interaction strength is increased, with 
the canonical signatures of such a transition in finite sys- 
tems. The functional r\ is depicted in Fig. [21 Each curve 
represents averaging over 10000 different disorder realiza- 
tions. In all cases 77 increases (the distribution becomes 
more "Poisson-like") as the interaction strength U is in- 
creased. This increase in rj becomes more pronounced the 
higher the filling factor, and tends to saturate at higher 
interaction strengths. At filling factors v > 1/4 (and 
v 7^ 1/2) our data suggest that 77 collapses onto a single 
scaling function up to ^-specific values where saturation 
takes place. 

Comparison with a SCHF scheme. It is known [f| 
that in zero-dimensional systems the low-lying many- 
body states are made of combinations of a finite number 
of Slater determinants (each made up of single-particle 
SCHF states). This motivates the use of such a basis 
to describe the exact |0) and |1) states. We invoke the 
participation ratio in the Hilbert space spanned by all 
possible Slater determinants: 
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Koopmans' variant of the HF approximation asserts that 
the effective single-particle states remain unmodified un- 
der a change of the occupancy of these states. As 
long as the HF-Koopmans-like picture holds 4] P = 1, 
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FIG. 1: The first excitation distribution -Pi(s) for different 
values of the interaction strength U. The Wigner-Dyson and 
Poisson distributions are plotted for comparison. 



FIG. 2: 771 (cf. Eq0J as function of the scaling variable Uv 
for various values of the filling factor v. Other than at low 
density [y < 1/4) and at half-filling the data seem to collapse 
on the same curve with ^-specific saturation values. 



whereas if this picture breaks down completely we ap- 
proach P ~ (jy ) , signaling derealization in Hilbert 
space @- 

Our results are depicted in Fig. |3J There we plot the 
disorder-averaged inverse participation ratio as function 
of the interaction for the ground-state and the first ex- 
cited state, with N = 4 and N = 7. It can be seen that 
1/(P) ex U 2 for weak interactions for both states. 

We note that there is some nontrivial physics taking 
place: the cusp of the curve 1/(P) vs. U for the first 
excited state, as it departs from the cx U 2 behavior (at 
U = 0.5 for the N = 4 case). This cusp becomes more 
pronounced as N becomes larger; the N = 7 case is a 
manifestation of the system becoming more localized ( in 
Hilbert space) with interaction strength exceeding U = 2. 
Studying the distribution of P (cf. Fig |3J reveals 
that as the interaction strength is increased, the original 
sharp peak at 1 is initially smeared, but the function then 
evolves into a doubly peaked distribution. 

In this context it is also useful to consider the quantity 
K = (l|cj v+1 CAr|0). We note that in the Koopmans limit 
K = 1, implying that both the ground and the first- 
excited states are given respectively by a single Slater 
determinant (they differ from each other by a particle- 
hole excitation). In that limit, as discussed above, the 
inverse participation ratio l/(-P) = 1. As we turn on 
the interaction 1 / (P) increases concomitantly with the 
decrease in K. The behavior of the average K is depicted 
in Fig El 

A two-level model. To further characterize the phe- 



nomena discussed here, we resort to the philosophy be- 
hind the Wigner surmise, and try to study the manifes- 
tation of the crossover in a truncated 2x2 Hilbert space. 
The latter is spanned by |0) and |1) (in a certain realiza- 
tion) . 

If the Koopmans-like scenario holds 0, 0, K = 1 
while (l|cjyCjv|0) = 0. Under such a scenario the Wigner 
distribution should hold: consider an arbitrary perturba- 
tion F = Y^ij a i 



(0|J10) 

<1|*|1> 
(1|F|0) = a N , N+1 , 



KN+l 

a ii + OCff+l,N+l 

i<N 



(6) 



corresponds (up to a constant matrix X x Yl i<N ctu sub- 
tracted, where X is the unit matrix) the 2x2 matrix: 



* H F — K oopmans 



OtN+l,N+l 



(7) 



The distribution of the HF states over disorder realiza- 
tions obeys RMT. Likewise, the distribution of {otij} will 
follow RMT (Gaussian) , and will lead to the energy spac- 
ing Ei — Eq obeying the WD distribution. 

Going now beyond the Koopmans limit, and employing 
the basis of the exact ground and first-excited states, the 



4 



0.15 




FIG. 3: The inverse averaged participation ratio 1/(P) and 
1/K = 1/(1|c] v+1 cjv |0), as a function of the interaction 
strength U for the ground state and the first excited state 
for the N = 4 case. The curves correspond to the best fit to 
1 + aU 2 , with a = 0.015 for the ground state and a — 0.22 for 
the first excited one. Inset: The same as the above for N = 7 
and a wider range of interaction of interaction strength. Here 
a = 0.034 for the ground state and a = 0.22 for the first 
excited one. 
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FIG. 5: The fluctuations in the matrix element, F a 
(a\F\(3), defined as {S 2 F a>0 ) = (F 2 ^) - (F a , 3 ) 2 , for the diag- 
onal (a=0) and non-diagonal (a ^ (3) elements as function of 
the interaction strength U. For TV = 4, 7 r s w 0.56U, 0.42C 
respectively. 



2-level matrix is 



F 



[1\F\1] (1\F\0) 
(1|F|0) (0|F|0) 



(8) 
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FIG. 4: The distribution function of the participation ratio, 
P, for the first excited state. 



The WD — > Poisson crossover is now manifest through 
F. We first subtract a constant, proportional to the unit 
matrix, to render the resulting matrix 6F statistically 
traceless. We are left to compare the fluctuations of the 
diagonal matrix elements, (5 2 F a ^) = (F 2 „) — (F Q ^) 2 
with a = j3, with those of the off-diagonal entries, 
(5 2 F a ^) with a ^ j3. (We note that the ensemble aver- 
age of the off-diagonal entries vanishes). This is shown in 
Fig. ISJfor two values of electron number, N. It is clearly 
seen, especially for the larger N, that as the interaction 
strength is increased (beyond U « 2 which is equivalent 
to r s = U /V^JAttv), the diagonal entries win over the 
off-diagonal, marking the suppression of level-repulsion 
and the cross-over to Poisson statistics. This cross-over 
appears to be quite sharp for our larger system. 
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